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ABSTRACT 

A galactic disk in a spiral galaxy is generally believed to be embedded in an extended dark matter 
halo, whicli dominates its dynamics in the outer parts. However, the shape of the halo is not clearly 
understood. Here we show that the dark matter halo in the Milky Way Galaxy is prolate in shape. 
Further, it is increasingly more prolate at larger radii, with the vertical-to-planar axis ratio monoton- 
ically increasing to 2.0 at 24 kpc. This is obtained by modeling the observed steeply flaring atomic 
hydrogen gas layer in the outer Galactic disk, where the gas is supported by pressure against the net 
gravitational field of the disk and the halo. The resulting prolate-shaped halo can explain several 
long-standing puzzles in galactic dynamics, for example, it permits long-lived warps thus explaining 
their ubiquitous nature. 

Subject headings: dark matter - Galaxy: Fundamental parameters - Galaxy: halo - Galaxy: ISM - 
Galaxy: kinematics and dynamics - Galaxy: structure 



1. INTRODUCTION 

Spiral galaxies are observed to have extended, nearly- 
flat rotation curves which indicate the existence of a dark 
matter halo. In fact, the radial distribution of the dark 
matter halo of a spiral galaxy is deduced from its ob- 
served rotation curve (Rubin 1983, Binney & Tremaine 
1987). However, the shape of the halo is not well- 
understood although it is expected to play a signifi- 
cant role in galaxy dynamics and evolution (Ryden 1990, 
Bekki & Freeman 2002). Different observational tracers 
such as tidal streams or thickness of the interstellar gas 
layer so far have given only the average halo shape, which 
in most cases, is found to be either oblate or spherical. 
In general, this is also the shape assumed in theoretical 
studies for simplicity. 

The vertical thickness of the interstellar gas provides 
an additional constraint, since it allows one to trace the 
force normal to the plane and hence the shape of the dark 
matter halo on the galactic scale. But this approach has 
been applied only to a handful of cases so far (Becquaert 
& Combes 1997, Oiling & Merrifield 2001, Narayan et 
al. 2005, Banerjee & Jog 2008, Banerjee et al. 2010, 
O'Brien et al. 2010). The vertical thickness of the HI 
gas is observed to increase sharply with radius in the 
outer Galaxy, which cannot be explained by a dark mat- 
ter spheroid of constant shape (Narayan et al. 2005), 
instead this requires a halo shape varying with radius 
(Kalberla et al. 2007). 

In this paper, we model the observed steeply flaring 
HI gas to constrain the shape of the halo, and show that 
the halo in the outer Galaxy is prolate. Further, the axis 
ratio varies with radius such that the halo is progressively 
more prolate in the outer parts with a maximum vertical- 
to-planar axis ratio of 2.0 at around 8 disk scale lengths. 
In Section 2 we describe the formulation of the problem, 
the results are presented in Section 3. Section 4 contains 
a discussion of related points, and Section 5 summarizes 
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our conclusions. 

2. FORMULATION OF PROBLEM 
2.1. Model for Disk Vertical Structure 

We have employed the multi-component model of the 
stars and gas of the galactic disk, in the field of the dark 
matter halo developed by Narayan et al. (2005). For sim- 
plicity, each disk component is assumed to be isothermal. 

The Poisson equation for an axisymmetric galactic 
system in terms of the galactic cylindrical co-ordinates 
(i?, (f), z) is given by 
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where pi with i = 1 to 3 denotes the mass density for 
each disk component, namely stars, HI and iJ2, and ph 
denotes the same for the halo. Here ^totai denotes the 
net potential due to the disk and the halo. For a flat 
or a gently-falling rotation curve, the radial term can be 
neglected as its contribution to the determination of the 
HI scale height is less than 10 percent as noted by earlier 
calculations (Narayan et al. 2005). 

The equation for hydrostatic equilibrium in the z di- 
rection is given by (Rohlfs 1977) 
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where < {v1)i > is the mean square random velocity 
along the z direction for the component i. We further 
assume each component to be isothermal for simplicity, 
so that the velocity dispersion is constant with z. 

Eliminating ^totai between eq. (1) and eq. (2), and 
assuming an isothermal case, we get 
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which represents a set of three coupled, second-order dif- 
ferential equations, one for each component of the disk. 
These are solved together to obtain the vertical density 
distribution of each disk component. This problem is 
solved in an iterative fashion, as an initial value problem, 
using the fourth order, Rimge-Kutta method of integra- 
tion. The two boundary conditions at the mid-plane i.e 
z = for each component are: 

Pi = {poh ^=0 (4) 

However, the modified mid-plane density (po)i for each 
component is not known a priori. Instead the net surface 
density Si(iZ), given by twice the area under curve of 
pi{z) versus z, is used as the second boundary condition, 
since this is known from observations. The HWHM (half 
width at half maximum) of the resulting model vertical 
distribution is used to define the vertical scale height. 

Our earlier study of the Galaxy based on this approach 
(Narayan et al. 2005), where a constant shape of the 
halo was used, showed that the observed flaring HI scale 
height distribution in the outer Galaxy can be explained 
by a spherical halo with the density falling faster than 
an isothermal case. This model, however, gave a total 
mass at the lower end of mass-range obtained by other 
techniques such as motions of satellites (Sackett 1999), 
which therefore questions its validity. 

2.2. Construction of the Dark Matter Halo Profile 

Here we propose and try an alternative viable idea, 
namely, a prolate halo with a shape varying with ra- 
dius. A prolate spheroid by definition has the vertical- 
to-planar axis ratio greater than 1, and it can be thought 
of as obtained by rotating an ellipse about its major axis. 
While a prolate halo fares better than a spherical halo of 
the same mass in explaining the flaring HI distribution, 
a single shape cannot explain the scale height data in the 
outer Galaxy over a large range of radii R = 9-24 kpc 
(Narayan et al. 2005, Kalberla et al. 2007). A halo that 
is progressively more prolate with radius is indicated to 
explain the observed steep flaring. 

At any R, an isodensity contour for a spheroid is de- 
fined by a constant "m" where = R'^ + z''^ jq^^ ; 
R! and z' are the co-ordinates of the points on the con- 
tour and denotes the vertical-to-planar axis ratio of 
the spheroid at R. Here R is the intercept on the mid- 
plane, that is, at z = 0. The mass of a shell bounded by 
two surfaces m and m -|- dm is a constant independent of 
the shape of the spheroid q^. ( see eq. 2.74 of Binney & 
Tremaine 1987 ), and this gives the following condition: 

QRPRiqR) = pB.{q = 1) (5) 

where pniqu) is the density along a prolate isodensity 
contour through R, and PR{q = 1) is the density along 
the corresponding spherical contour with radius R. Thus 
if the original spherical shell were instead taken to be a 
prolate spheroidal shell by construction, its density will 
be lower by a factor qn. Thus the vertical force near the 
mid-plane, a crucial determinant of the HI scale height, 
will be lower and this could qualitatively explain the HI 
flaring that is observed. 



2.3. Solution of Equations for a Prolate Halo 

We start with the Galactic mass model obtained by 

fitting the observed rotation curve (Mera et al. 1998) 
which gives the pseudo-isothermal spherical halo density 
distribution at a radius r to be: 



where po = 0.035 M© pc~^ and rc = 5 kpc are respec- 
tively the best-fit central density and the core radius for 
the halo. This determines the mass within a spherical 
shell at each r — R. If the shape were taken as prolate 
instead, the density will change as per eq. (5) while the 
mass in the shell remains constant. One can ensure that 
the fit to the rotation curve still remains valid with devia- 
tion of only a few percent since the rotation velocity has a 
weak dependence on the shape of the halo mass distribu- 
tion (Sackett & Sparke 1990). In contrast, a small change 
in the halo shape of a given mass can have a striking ef- 
fect on the vertical scale height. We exploit this idea to 
constrain the variation in the halo shape, q^ , with the 
radius R. In a sense, each shell acts independently. We 
model qn as a simple second-order polynomial in R i.e 
qR = 1. + ai{R —9) + a2{R — 9)^ over a range R = 9 - 
24 kpc. Here Q!i,a2 are free parameters. Combining this 
trial qR with eq. (5) and (6) gives the halo density, ph, 
which is used as an input parameter in solving eq. (3). 

The main aim of this paper is to investigate whether 
a prolate halo explains the observed HI flaring within 
the purview of the above mass model. We next calculate 
the HI scale height versus radius theoretically by solving 
eq. (3) numerically with the procedure as outlined in 
Section 2.1. The input parameters such as the stellar 
surface density, and stellar and gas velocity dispersions 
are taken to be the same as in our earlier work on the 
Galaxy (Narayan et al. 2005). The surface densities of 
HI and H2 are taken from observations (Wouterloot et al. 
1990). The resulting HI scale height values obtained by 
solving eq. (3) for different trial are qn values are flt to 
the observed data (Wouterloot et al.l990) over the range 
R = 9-24 kpc. The halo is taken to be spherical at i? = 
9 kpc. Since the error bars on the data are not available 
(Wouterloot et al. 1990), therefore we assume error bars 
of 5% and obtain the best-flt values of the parameters 
ai,a2 using the method of least squares. These are: at = 
0.020 and a2 = 0.003 and the range ai = 0.010 - 0.080 
and Q!2 = 0.001 - 0.004 gives reasonably good flt as well 
(Bevington 1969). 

3. RESULTS 

The plot of the resulting vertical scale height versus i?, 
and the comparison with the observational data which 
were corrected for the Galactic warp (Wouterloot et al. 
1990), is shown in Figure 1. The corresponding best- fit 
qn increases from a value of 1 corresponding to a spher- 
ical halo at i? = 9 kpc to a prolate halo with qij = 2.0 
at ii = 24 kpc. We stress that the qn value denoting the 
halo shape as obtained here is a local property, and our 
approach has allowed us to obtain a radial variation in 
it. In contrast, studies involving tidal streams (Ibata et 
al. 2001, Law & Majewski 2010) and other observational 
tracers assumed the halo shape to be constant, and found 
it to be typically spherical or oblate {qR < 1) (Sackett 
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1999). There arc a few exceptions involving studies of mi- 
crolensing (Holder & Widrow 1996) and stellar streams 
(Helmi 2004) which yield a prolate halo, but again with 
a global, constant shape qr < 1.5 within i? = 60 kpc. 

Figure 2 gives the isodensity contours (on the R-z 
plane) denoting this increasingly prolate halo in the outer 
Galaxy, and the inset shows the resulting qn versus R. 
We emphasize that our study traces the halo shape vari- 
ation over a large radial distance upto 24 kpc, that is 8 
disk scale lengths (with the disk scale length being 3.2 
kpc as in the Mera et al. (1998) model) or twice the 
optical radius of the Galaxy. Surprisingly, a prolate halo 
with a maximum qn — 2.0 at i? = 24 kpc is sufficient 
to explain the steep rise in the gas thickness by a fac- 
tor of 6 seen over a radial range of 9-24 kpc. This is 
because the density in a prolate shell at a given R is in- 
versely proportional to qa (see eq. 5), while the vertical 
force due to gravity and hence the resulting scale height 
depends in a nonlinear way on the density near thc^ mid- 
plane (Banerjee & Jog 2007). Thus a small decrease in 
the midplane density is sufRcient to explain the observed 
sharp increase in scale heights. 

4. DISCUSSION 

1. Dependence on gas velocity dispersion: The 
value of gas dispersion plays a crucial role in these cal- 
culations, but it is not easy to measure and depends on 
the fraction of HI in the different phases. In the outer 
Galaxy the ratio of the 21 cm emission and absorption is 
shown to be remarkably constant upto R=25 kpc which 
implies a constant ratio of the warm and cold phases 
of HI (Dickey et al. 2009). However, the pressure sup- 
port for the gas is mainly from non-thermal or turbulent 
motions which dominate the thermal velocity dispersion, 
and Dickey et al. (2009) do not give the non-thermal 
velocitites. Hence here we use the measured dispersion 
at the solar location of 9 kms~^ (Malhotra 1995). This 
is expected to decrease at larger radii and saturate to 
7 kms^^ (see Narayan et al. 2005). These values were 
used as input parameters in our calculations. These val- 
ues are in good agreement with the radial variation seen 
in external galaxies (Kamphuis 1993, Tamburro et al. 
2009). 

2. Mass of the prolate dark matter halo: In this 
paper we start with the screened isothermal mass model 
of the Galaxy as given by Mera et al. (1998), which 

has a mass of 10^^ Mq within R < 100 kpc. Due to the 
construction adopted in this paper, the mass in each shell 
is conserved as an initially spherical shell is stretched into 
a prolate shell (eq. [5]). Thus the total mass of the halo 
is conserved. Hence the total mass of the prolate halo in 
our model is the same as in the Mera et al. mass model. 
We note that in our model we have only dealt with the 
region of i? < 24 kpc. 

This total mass is in a good agreement with the value 
obtained using different observational tracers such as mo- 
tions of satellites, and high velocity stars (Sackett 1999, 
also see Gnedin et al. 2010, McMillan 2011 for more re- 
cent estimates). Thus the current model overcomes the 
problem with the earlier model (Narayan et al. 2005) 
which gave a total halo mass within 100 kpc that was 
three times smaller than the above value. Recall that 
this was the motivation for trying the variation in the 



halo shape with radius (Section 2). Thus our model ex- 
plains the observed HI scale heights and also gives a total 
mass in agreement with typical values in the literature. 
3. Asymmetry in gas scale heights: The scale 
height distribution is known to be asymmetric in the two 
halves of the Galaxy (Levine et al. 2006, Kalberla et al. 
2007) being higher by a factor of about two in the north- 
ern galactic hemisphere than in the south. We have only 
modeled the northern data (Wouterloot et al. 1990) here 
for simplicity to check if our proposed model of a prolate 
dark matter halo with a radially-dependent shape could 
explain the steeply increasing HI scale height data. To 
model the southern region, we would need the observed 
HI surface density values as input parameters. Levine et 
al. (2006) measure the HI scale height for both north and 
south, but only the net value of the gas surface density, 
averaged over the north and south is given. 

Recently, Kalberla & Kerp (2010) have given observed 
HI surface density values separately for the north and 
south. Applying our model to the southern data alone, 
we find that the best-fit qa values are significantly lower. 
A prolate halo with a radially-increasing axis ratio is 
still preferred over an oblate or spherical case, but with 
smaller values of ai = 0.009 and 0:2 = 0.0 (compared to 
0.02 and 0.003 respectively for the northern data. Section 
2.3). These give; a maximum axis ratio of 1.14 at i?=24 
kpc as compared to qa = 2.0 obtained for the north. 
This low value of qR is sufficient to explain the data in 
the south, since the scale heights are lower and the rate of 
fiaring beyond 18 kpc is much less steep in the south than 
in the north. Thus the dark matter halo of the Galaxy 
has a prolate shape but its value is different in the two 
hemispheres, being more prolate in the north. A more 
realistic future treatment should consider simultaneous 
modeling of north and south HI scale height data. 

There is a limitation in applying this model farther 
out in the plane since the gas dispersion values are not 
well-known. With this caveat in mind, if we model the 
data the range 9-40 kpc (Kalberla & Kerp 2010) by 
extrapolating the above input parameters, the best fit 
gives a maximum qji = 4.5 at i? = 40 kpc. We have 
not modeled the inner Galaxy (i? < 9 kpc) here, since 
the halo contribution is shown to be small in this region 
(Narayan & Jog 2002), hence the scale heights cannot be 
used to constrain the halo shape. 

4- Comparison with shapes from cosmological 
simulations: Interestingly, the prolate-shaped dark 
matter halo obtained here by modeling the HI vertical 
scale height data agrees well with the general trends seen 
in cosmological simulations. The latter give a range of 
halo shapes with a preference for a prolate shape (Bailin 
& Steinmetz 2005, Bett et al. 2007) but with a lower 
vertical-to-planar axis ratio. However, these are mea- 
sured at scales of about 100 kpc (Bailin & Steinmetz 
2005) which are much larger than that of a galactic disk, 
hence the quantitative comparison is not meaningful. 
The predictions of our model can be checked and con- 
firmed by the signatures of other local tracers in future 
observational studies such as GAIA. 
5. Dynamical implications of a prolate halo: A 
prolate-shaped dark matter halo as obtained here has im- 
portant implications for many classic problems in galac- 
tic dynamics. For example, a prolate halo would cause a 
lower differential precession and hence can support long- 
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lasting warps (Ideta et al. 2000), this can explain why 
warps are commonly seen. Second, a prolate halo can 
naturally explain (Helmi 2004) the fact that the satel- 
lites are seen to be limited to the polar plane normal to 
the galactic plane, as observed in many galaxies, known 
as the Holmberg eflFect. The dynamical implications of a 
prolate halo, in particiilar one which is increasingly more 
prolate with radius, deserve a detailed study. 



5. CONCLUSIONS 

By modeling the flaring HI gas distribution, we have 

shown that the shape of the dark matter halo is prolate 
in the outer Galaxy over a radial range of i? = 9 — 24 kpc, 
where it can have possible observable dynamical conse- 
quences. We treat the halo as a set of spheroidal shells 



which have progressively more prolate shapes with in- 
creasing radii. Thus, wc obtain the local shape of the 
halo and show that the maximum vertical-to-planar axis 
ratio is 2.0 at i? = 24 kpc. These results are in contrast 
with most of the earlier work involving various obser- 
vational tracers, which gave either a spherical or a flat- 
tened, oblate halo with a constant shape. However, our 
results agree with the trends from cosmological simula- 
tions which tend to favor a prolate halo. A prolate halo, 
in particular one which is increasingly more prolate at 
larger radii, has important implications for galaxy dy- 
namics and evolution, which need to be studied further. 
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Fig. 1. — The calculated vertical scale height for the atomic hydrogen gas, HI, (solid line) and the observational values (Wouterloot et al. 
1990) (filled circles) versus Galactocentric radius R. The theoretical curve is the best-fit case, and corresponds to a dark matter halo which 
is progressively more prolate with radius. In the radial range studied, the halo is found to be most prolate with the vertical-to-planar axis 
ratio, qji = 2.0 at ii = 24 kpc. 
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Fig. 2. — The resulting best-fit prolate-shaped isodensity contours of the dark matter halo on the R-z plane, with vertical-to-planar axis 
ratio increasing with radius as depicted in the inset. This clearly brings out the increasingly prolate geometry of the halo shape. The 
densities of the successive contours moving radially outwards are 8.3, 5.2, 3.5, 2.5, 1.9, and 1.5 in units of lO""^ Mq pc""^. 



